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ABSTRACT 

Large single crystals of benzene (C6H6) were aligned by an 
optical method. The aligned crystals were pulsed with ultrasonic 
waves at a frequency of approximately 10 megacycles per second using 
quartz crystal transducers. The velocity of the ultrasonic waves 
was measured by comparing the echo transit time with that of simul- 
taneously generated waves moving through a variable length path in 
a medium in which the velocity of sound was known precisely. 

Velocities of longitudinal and transverse waves along each of 
the principal axes of the orthorhombic crystal were used to compute 
the main diagonal terms of the “4 ጎ matrix. The off-diagonal terms 
were then determined by velocity measurements along arbitrarily 
selected directions. 

Velocity measurements were made and the elastic constants 
computed over the temperature range 170°K to 250%. 

Absorption in the crystal was measured over the same temperature 
range for each of the directions used for velocity measurements. The 
absorption was determined to be aelotropic. 

The writers wish to express their appreciation for the assistance 
and encouragement given them by Professor ሀ. B. Wilson of the U. 5. 


Naval Postgraduate School in this work. 
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INTRODUCTION 

Benzene, the fundamental aromatic hydrocarbon, possesses an 
orthorhombic crystalline structure in the solid state. 

The low melting point of benzene (5.2°C) renders the use of 
standard diffraction methods for the determination of orientation 
difficult. Benzene crystals are, however, transparent to visible 
light and possess a high positive birefringence; thus their optical 
properties provide a suitable means of alignment using simple optical 
methods. 

The formulation of the strain energy equation for an orthorhombic 
crystal results in nine independent elastic constants(1). By measure- 
ment of the velocities of longitudinal and transverse ultrasonic waves 
along the three principal axes, the six main diagonal terms of the Ci; 
matrix can be computed. The remaining three elastic constants must 
then be determined by velocity measurements along arbitrarily selected 
directions, using the previously determined values for the main 
diagonal terms. 

Using the techniques for alignment, pulsing and velocity measure- 
ment described in this thesis, the velocity data was obtained for the 
computation of the elastic constants over the temperature range 170°K 
to 2500K. 

In obtaining the velocity data, the opportunity was taken to 
observe and record the absorption of the ultrasonic pulses at a fre- 
quency of 10 megacycles per second as a function of crystal orientation 


and of temperature. 





CRYSTAL GROWING AND CRYSTALLINE STRUCTURE OF BENZENE 





The large single crystals, necessary in order to follow the align- 
ment and cutting procedures developed, were grown in a melt of reagent 
grade benzene. The melt was contained in a sealed container which was 
maintained at a temperature of approximately 5°C. After degassing, the 
bath was left at this temperature for approximately two weeks during 
which time the crystals formed. 

The conglomerate mass of crystals was then broken apart and large 
single crystals of irregular shape were obtained and stored at a 
temperature of about -10°C until required. 

Crystalline benzene is of the orthorhombic class, space group 


Pbca, the unit cell having the dimensions (2): 


a= 7.20219 
b = 9.1471 AO 
Cea Oe eno ሚመ ር 


Optically the crystal is transparent with indices of refraction (3): 
ጋ 11 
D = = 0.00) 
ia an? 
The crystals were thus very suitable for examination by optical 
methods. 
CRYSTAL ALIGNMENT 
In order to obtain the necessary data to solve the strain energy 
equation and thereby compute the elastic constants, crystals were 
required which were accurately aligned along each of the principal 
axes together with at least one sample which could be aligned along 


2 





some arbitrary axis. 

The method of alignment used consisted of examining the raw ir- 
regularly shaped crystals under a conical beam of polarized light. 

For this purpose, a conoscope was constructed in the form shown 
schematically (Fig. 1). For actual examination the instrument was 
located in a refrigerated compartment maintained at a temperature of 
about -10°C. 

In the raw state, since there seems to be no apparent preferred 
direction of growth for benzene crystals, initial examination was made 
with the crystal immersed in a glass bottomed tank containing a mixture 
of benzene (CgH¢) and methyl salicylate (CgHgO 3) - This mixture was in 
the proportions of lO parts methyl salicylate to 300 parts benzene by 
volume. The refractive index of this mixture is approximately 1.52, 
close to that of benzene at ~10°C. The actual mixture proportions were 
found to be not critical. Changes occur due to evaporation and the 
solubility of the crystalline benzene in the mixture. Fortunately 
the mixture is close to the solubility limit of benzene in methyl 
salicylate, so that loss of crystal by dissolution is small. 

The tank was positioned in the conoscope so that the bottom of the 
fluid was at the focal point of the convergent lens, with the divergent 
beam spreading up through the fluid. 

In this beam the raw crystal was examined for the location of the 
optic axes, the bisectrix of which locates the b axis of the crystal 
Piz. 8). 

Examination was carried out by rotating the crystal, held in a 


pair of tongs, in the beam to determine two extinction planes. Since 
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the normals to these extinction planes form an orthogonal system with 
one of the optic axes, this axis was visible at 90 degrees to both these 
normals. 

This alignment was retained by transferring the crystal to a special 
clamp (Fig. 2), designed to accommodate a range of shapes and sizes of 
crystals and at the same time permit a more accurate alignment. The 
clamp consisted of a short length (about 5cm) of 10cm diameter aluminum 
tubing into which were mounted in gimbal fashion two aluminum rings. On 
the inside of the inner gimbal were secured four strips of spring steel 
bent into the center forming a spring collar which retains the crystal. 
Spring pressure was just sufficient to hold the crystal firmly in 
position. 

Mounted in this clamp, the crystal was returned to the bath with 
the base of the clamp forming a reference face resting on the glass 
bottom of the bath and was thus in a plane normal to the optic axis of 
the conoscope. 

With the freedom afforded by the gimbals, the crystal could be 
rotated with respect to the optic axis of the conoscope to bring the 
bisectrix into coincidence with this axis. Since the crystal was held 
firmly in the clamp, the interference pattern (Fig. 3) could be clearly 
observed. 

The alignment was preserved by removing the clamp and crystal from 
the bath and, using the base of the clamp as reference, melting paral- 
lel faces on the crystal, which were finally covered with glass slides. 

A check on this alignment was made by mounting the crystal on a 


coniometer attachment, which could be fitted to the conoscope. 
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Small corrections to the alignment could be made by melting new 
faces on the crystal, the corrections being computed from the goni- 
ometer angles using the knowledge of the refractive indices of benzene. 

Alignment to the other principal axes was accomplished by rotating 
the goniometer 90 degrees in the plane of the optic axes (190) eram a 
plane normal to the plane of the optic axes (001). After rotation, new 
parallel faces are melted ja die crystal and glass slides secured to 
preserve the particular orientation. 

Any other orientation could be obtained by appropriate rotation 
of the gonicmeter. In all cases, a knowledge of the direction of the 
optic plane was preserved by appropriately marking the glass covering 
slides. 

MOUNTING OF CRYSTALS 

For the purpose of ultrasonic pulsing, it was necessary to mount 
the crystals in such a manner as to permit electrical connection of a 
radio frequency source and a receiver to a transducer, which was in 
mechanical contact with the crystal. At the same time, it was 
necessary to ensure that the crystal orientation was maintained and 
the effects of sublimation eliminated. The crystal holder (Fig. h) 
was therefore designed and manufactured to meet these requirements. 

To prepare a crystal for mounting in this holder, it was trans- 
ferred to a glass plate and frozen in place using the preserved face as 
reference; the desired axis was thus maintained normal to the glass 
plate. The cylindrical shape required to fit the holder was then 
obtained by melting the periphery of the crystal, using metal cylinders 


maintained at room temperature to melt the crystal. When the desired 
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cylindrical shape was obtained, the holder was slipped over the sample 
and frozen in place. During this process, the sample and holder were 
maintained at refrigerator temperature. 

With the crystal in the holder, the end faces were finished 
parallel to the holder faces by polishing on fine emery cloth. The 
end faces were then sealed with 1/) mil clear Mylar for the free end 
and 1/ mil aluminum coated Mylar for the end receiving the trans- 
ducer, using glycerine as a sealing agent. The end seals were secured 
in place using the holder end plates. 

The transducer was mounted on the aluminum coated Mylar using 
Nonaq stopcock grease. These bonding materials were found satisfactory 
for the transmission of both shear waves and longitudinal waves over 
the temperature range for which data was obtained. 

The transducers employed were 0.5 inch diameter quartz crystal. 

X cut crystals were used for the generation of longitudinal waves; 
Y cut for shear or transverse waves. 
PULSING OF CRYSTALS 

For pulsing, the crystal holder was mounted in the pulsing 
assembly (Fig. 5). This assembly provided a secure support for the 
holder and at the same time permitted connection of the transducer to 
the radio frequency source. The ground connection was made through 
the body of the assembly and the crystal holder, the other connection 
being made through the spring loaded electrode. The method of mounting 
the holder allowed adjustment of the electrode pressure on the trans- 
ducers a light pressure was found best for the generation of longitudi- 


nal waves and a heavier pressure for transverse waves. 
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Radio frecuency pulses were generated by an Arenberg pulsed 
oscillator, Model PG 650-C, connected by low impedance coaxial cable 
to the pulsing assembly. The return echoes after reflection from the 
free end of the crystal were amplified by an Arenberg wideband ampli- 
fier WA-600 and displayed on a Tektronix Type RML5A oscilloscope. A 
Tektronix Type RM 181 time mark generator was used for control of pulse 
repetition rate and for external sweep control of the oscilloscope. 
Figure 6 shows a schematic arrangement of the pulsing system. 
MEASUREMENT OF VELOCITY 

To obtein accurate measurements of the velocity of the ultrasonic 
waves through the crystal, a modification of the method of velocity 
measurement described in (h) was used. Figure 7 shows the construction 
of the velocity comparator designed and built for this modified method. 
In the velocity comparator, pulses generated by the quartz crystal 
transducer immersed in the liquid medium were reflected off the rigid 
boundary and received at the transducer. The path length and thus the 
transit time of the pulse was made variable by altering the position 
of the boundary by means of the micrometer screw. 

In making a measurement, the pulse trains from both the comparator 
and the crystal were superimposed on the oscilloscope and the path 
length adjusted to bring corresponding echoes into coincidence. Pre- 
cise matching could be obtained by observing the interference effects 
as the echoes were brought into coincidence, since the comparator and 
benzene transducers were pulsed in parallel. 

The velocity in the crystal was then computed from the known 


velocity in the liquid medium by the relation: 
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where: 1, = initial micrometer reading corresponding to zero 
path length. 


1, = length of path through crystal. 

1 = micrometer reading at coincidence. 
Vu = sound velocity in the liquid medium. 
Ve = sound velocity in the crystal. 

The liquid medium used was a mixture of 100 parts distilled water 
to 22 parts absolute ethyl alcohol by volume. This mixture transmits 
sound waves at a velocity given by Willard (5) which is of a magnitude 
suitable for comparison with the velocities of benzene crystals and is 
substantially independent of temperature in the region of 25°C. 
Temperature control of the device was, therefore, unnecessary. The 
temperature, however, was monitored and corrections made for the 
small variations which occurred. 

Corrections to the measured length of the crystal were necessary 
to take into account the effect of thermal contraction of benzene. 
These corrections were made from the data on the thermal expansion of 
benzene given by Cox (2), and by Euchen and Lindenberg (6), and applied 
as a percentage decrease in the length of the crystal. 

The accuracy of the velocity measurements was dependent on two 
dominant factors: 

(a) Measurement of crystal length: this length was measured to 
«001 cm and, after correction for thermal contraction, was considered 


accurate to within .01 cm for the 3 cm path length. 
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(b) Pulse matching between the crystal and velocity comparator; 
the two 10 mc pulses were matched to within one half of one wave length, 
resulting in a possible error of ‘05 sec in 12 u sec. 

The possible errors introduced by other factors were not of 
significant magnitude in comparison to these two factors. Based on this 
estimate, the velocity error was determined to be less than 1%. 

Measured velocities, using the same crystal, were reproduced to 
within 1% over the temperature range 250°K to 170°K. 

TEMPERATURE CONTROL AND MEASUREMENT 

Velocity measurements were made over the temperature range 250°K 
to 170°K by cooling the crystal mounted in the pulsing assembly using 
liquid air. In order to maintain an even cooling rate the sealed 
pulsing assembly was inserted in a silvered Dewar containing liquid 
air. Contact with the liquid air was made through a heat leak, which 
permitted control of the cooling rate. 

Temperature was measured using a Copper-Constantan thermocouple 
placed in contact with the crystal holder. 

TEMPERATURE RANGE 

The upper limit 250°K was predetermined by the need to preserve 
the crystal from sublimation and also ensure the maintenance of a 
satisfactory bond between the transducer and the crystal. 

Bond failure occurred in the temperature range 150°K to 170%, 
thus determining the lower limit of 170°K. The precise nature of the 
failure was not ascertained but it was believed that differential con- 
traction of the crystal and the Mylar end seal resulted in separation at 


the glycerine bond producing a void between the transducer and the crystal. 
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DETERMINATION OF ELASTIC CONSTANTS 
For a crystal of orthorhombic symmetry, Group 3, the nine inde- 


pendent elastic constants may be arrayed in accordance with Cady (1), 


as shown: 
“ገዝ Cj2 ‹ገን 0 ዐ ዐ 
C12 Cop Co3 0 ዐ ዐ 
GE C53 C33 ዐ O O 
O O O Du O 0 
O 0 0 ዐ “55. 0 
ዐ 6 O O O C66 


Since it was possible to align and mount crystals oriented along 
each of the principal axes, the main diagonal terms Cs; of the matrix 
were computed directly from longitudinal and transverse wave velocity 
measurements along each of these axes. 

The crystal orientation, measured velocity and associated elastic 


constant determined by this method were as given in Table 1. 
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Note: For axial velocities, the first letter indicates the axis of 


propagation, second letter indicates axis of wave displacement. 
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To determine the off-diagonal terms Ci > tne transformation 
equations (App. I) were used. Crystal orientations were chosen so that 
each off-diagonal constant could be computed in terms of a measured 
velocity and certain of the previously determined main diagonal 
constants. For the orthorhombic crystal the desired form of the 
equations was obtained by considering a rotation in an axial plane to 
an angle of 5 degrees from two of the principal axes (Fig. 8). This 
produced three different orientations giving three independent sets of 
equations, each of which could be solved for one of the three off- 
diagonal terms. 

The three independent sets of equations were obtained by sub- 
stituting the appropriate values for the direction cosines 1, m, and 
n in the equations as demonstrated in Appendix I. 

Substitution of these particular values reduced two of the Ci; 
terms to zero. The remaining terms when substituted into the secular 
eauation for irrotational plane waves traveling in the X’ direction 
gave rise to the simplified determinant solution: 

(Ci - pve ) (C55 = pve ) - 6445 = 0. 
This may be reduced to: 

pesa let verd) o MES | 
where the C” terms are defined in terms of the elastic constants and 
the direction cosines (App. I). This results in an equation in terms 
of the longitudinal velocity and previously determined elastic constants 
(044) 


ES - É vers L(C: : የ DE y 2 (Cs; C ye C 
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The crystal orientation, measured velocity and associated elastic 


constants were as given in Table 2: 





Table 2 
Direction of ። .: - Velocity -ር.. - ር.. TS 
Wave Propagation Measured AAA _ JJ Kk | 
O Ye Ge Du Ca 066 
Si. Voll C23 2 Cuh 


All of the data pertinent to the computations and the actual 
values of the elastic constants have been presented in both graphical 
and tabular form, page 21-26. 

ERROR IN ELASTIC CONSTANTS 

The sources of error in the determination of the elastic constants 
from the measured velocities could be divided into two categories - the 
error in crystal alignment, and the error in computing the elastic 
constants from the velocity data. 

The error in crystal alignment depended not only on the accuracy of 
the optical alignment, but also on the method of preserving this align- 
ment during the mounting procedure. Since for all alignments other than 
that along the b axis, the alignment of the crystal could not be 
rechecked after mounting, no correction could be applied. 

Another source of alignment error was introduced in the alignment 
of the shear transducer to produce shear waves along the axes. To 
eliminate this error, only longitudinal wave velocities were used to 


determine the off-diagonal elastic constants. 
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For the crystals not aligned on a major axis, a further source of 
error was introduced by the anisotropic thermal contraction, which 
produced a continual minor change of alignment throughout the tempera- 
ture range. Because of other larger sources of error present in the 
determination of the off-diagonal constants, no correction was made 
for this error. 

The specific volume of benzene was considered to be 0.9125 cc/gm 
for the temperature range 250°K to 170°K, based on data summarized in 
Andrew and Eades (8). Possible errors in this data were considered 
negligible in comparison with the other errors present. 

Since it was not possible to measure any alignment errors present, 
an estimate of the accuracy was obtained by comparison of velocities 
which determine the same elastic constant. For the constants determined 
by longitudinal axial velocities, the error was estimated at less than 
2%, based on velocities obtained from different crystal samples. For 
the constants determined by shear velocities, the comparative 
velocities obtained also differed by no greater than 2%. 

The values of the diagonal elastic constants (C4 ;) determined from 
axial velocities were required in the solution of the off-diagonal 
elastic constants (Cj;j); thus, any error present in these values was 
unavoidably included in the off-diagonal values. This resulted ina 


possible overall error in the off-diagonal elastic constants of up to 






10%. 
Table of Estimated Error in Elastic Constants 
Elastic Constants E | “Velocities Used t 
C11, Co, 38 % ብ “Longitudinal, Axial | less than 2% 
Chi, Coo, Dee Shear, Axial approximately 2% 
ር lene 24 Longitudinal, Non-axial up to 10% 
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MSASUREMENT OF ABSORPTION 

Although measurement of absorption was not the primary objective, 
the method used to pulse crystals to obtain velocity data was readily 
adaptable to providing absorption data. 

Absorption measurements were, therefore, made in conjunction with 
velocity measurements by comparing the amplitude of successive pulse 
echoes. The comparison was made by matching the leading echo to the 


output of the time marker generator displayed on the second channel of 
the oscilloscope. The time marker generator output was fed through an 
attenuator which permitted an adjustment of amplitude in steps of one 
decibel. 

In measuring the attenuation of successive echoes, the time marker 
was matched to the leading echo and the attenuator then adjusted until 
the time marker amplitude was the same as that of the succeeding echo. 


The attenuation between echoes was, therefore, given directly as the 
attenuator setting in decibels. These measurements were repeated for 
as many echoes as appeared in a given pulse train. 

All absorption measurements were made at a frequency of 10 mega- 
cycles per second using 0.50 inch diameter quartz crystal transducers. 
The accuracy of the measurement. was extended to within one half of one 
decibel by logarithmic interpolation of the superimposed waveforms. 
DETERMINATION OF ABSORPTION 

The average of the attenuator settings was then converted to 
absorption in nepers per centimeter by a calculation based on the 


measured length of the crystal: 


Absorption (nepers/cm) = attenuation (25) x 8.68 nepers 


24 (240 
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The ultrasonic absorption computed from the measurements of 
pulse attenuation has been presented graphically on page 30 and in 
tabular form on page 31. 

A11 of this data was based on the attenuation between first and 
second echoes only, thus involving a total path length of six centi- 
meters in the crystal. For this length of path, the effects of 
diffraction were small and within the errors of the actual measure- 
ments, thus no correction was made for attenuation due to diffraction. 

The absorption computed varied from greater than 0.6 nepers per 
cm. to less than 0.1 nepers per cm. the variation being a definite 
function of crystal alignment in relation to the direction of ultra- 
sonic wave propagation. The absorption for all waves propagating 
along the c axis was significantly less than for waves propagating 
perpendicular to this direction. 

Some variation of absorption with ገሽ ce was observed for 
certain of the crystal directions. In general this variation showed 


a decrease in absorption with a decrease in temperature. 
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ULTRASONIC ABSORPTION 
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ULTRASONIC ABSORPTION IN BENZENE 


Direction Temperature Attenuation Absorption 
of Range nepers 
Propagation OK db/cm. per cm. Mode 

100 250 - 170 5.3 0.61 Shear (001) 

010 250 - 230 283 0.38 Shear (001) 
230 = 170 በ. 0.19 

110 250 - 230 640 0210 Long. 
250 SO 3.5 0.40 

101 240 - 225 TEN 0.54 Long. 
225 - 210 』 。2 0.48 
210 - 195 ር. 0.ከዐ 
195 jan 170 0 0.35 

010 250 = 230 በ 0.53 Long. 
230 - 210 e 0.18 
21.0 ro 200 Ba 0.13 
200 - 170 3.0 0.35 

100 250 - 2h0 383 0238 Long. 
240 = 220 3,0 0:35 
220 - 170 の 0524 

001 250 - 210 2.0 0.23 Long. 
2ከ0 = ጊ70 E 0.19 

011 250 = 215 1 0.20 Long. 
216 = 205 1.5 Deeg 
205 = 170 ee 0.15 

001 250 - 210 223 0.26 Shear (010) 
210 = 190 ims 0.15 
190 - 170 057 0.08 
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APPENDIX I 
The method develcped by J. R. Neighbours (7) for the determination 
of elastic constants is based on velocity measurements of irrotational 
plane waves propagated along a specified reference axis ESA where the 
reference axes are oriented so that the Y’ axis lies in the 001 (XY) 
plane of the crystal. 
Using the measured velocity data, the elastic constants are 


determined from the solution of the secular equation: 


( ュ -ZV2) 016 Gé 
Cag (C2 - pV?) Cee = 0 
015 556 (55 -ዖ"") 


where the primed constants are defined in terms of the crystal elastic 


constants by: 
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+ (mnC23 + {mna-2(a2-n°)0)}, = (mna~@(a@=n?)Coe -2/mna-2(( ?-me)Cgg 
mere a? = (E + m? 
De 





2 


/ 
and the 2 m, and n terms are the direction cosines relating the X 
direction to the crystallographic axes (a,b,c). 
For a wave propagating along a crystallographic axis, the direction 


cosines become one or zero, which results in: 


7 
ዐ4 4 
ሥ 


Cii = Ci; 


= 0 


So that the solution is in terms of 


pve = (411 

ሥ 
Example: for el Eeer q ዖ"ሜ« 
m=n = 0 Ger = C55 * ዖሻ"3እ 


Sa a, 

For a wave propagating in an axial plane at 5 degrees to two 
crystallographic axes, the direction cosines become zero or 12/2. 
For the 110 case, where the wave is propagating at 15 degrees to the 
a and b axes, and perpendicular to the c axis: 


-ር= በ ። の C14 = ፳(011 + Coo + Jee + 2C] 2) 


1 
O C66 = «(Cy + 622 >2012) 
a = 1 C55 = SO A Cgg) 
C16 = #(C22 -C11) 


e 


/ 
Cis = “ሮይ = 0 


thus the determinant becomes: 


ュー PV?) C16 0 
C16 (cee - Pp) 0 = 0 
0 O (0ፊ - pv?) 


or (Cee = pv?) | (ary - ዖ"")(ዐሬ - ዎ"") - 462] =0 


SE 





which has solutions: 


ሠ 


(a) „ve = ウロ = 


Ue 


LO 7 655) 


At 


Al 
ገ E / 3 £ ሥ = 
(5) pv? = C11 + 066) 4 Ken - 066)? + 0162 | 


The solution for the longitudinal velocity (V/) ís: 


L 
E 4 ei “ “ 2 
pv, © = 361, + 068) + (40077 - 066)? + 0162] 


wir 


= 401] #092 # 2066) +3 [(012 + 066)? + Ž(Cop D? 


Separating the off-diagonal elastic constant 


A JL 
8 ኃ 2 
Cio - 172% - 3(0;1 * 022 * 2ርሬፊ)]- - 2(022 - e?) - C66 


results in a solution for the off-diagonal constant in terms of the 
longitudinal velocity and the previously determined diagonal constants. 
The solutions for the other off-diagonal constants are analagous: 
for the 101 Case C3 = pr? - HC + C33 + 20e) - 4039-02)” - Ccg 
- 2 7 
for the 011 Case Co っ ka? - (Cpo + C33 + Ce) | - X(033-022) 2) à - Ch, 


፲ቱ may be noted that since the solution is of the form: 


(Ci; + Ch) = a known positive constant (k2), 
this allows for the possibility of two solutions: 
Ci; = k - Oy 
EE 


or Ci; 
The choice between solutions was made on the basis of the stability 
requirements for a crystal of orthorhombic symmetry. From Neighbours (9), 
this requires that the principal minors of the matrix must all be 


positive. The matrix was therefore expanded into terms of the principal 


3h 








minors, and the following four inequalities obtained: 
013022 › 612" 
611033 > 013° 
522033 >> 023" 
2:1 (022፡33-023") + 2032003013 > 0992033 + 043059 
These inequalities were then used to show that the solution giving a 


large negative value was not acceptable. 


BG 
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